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ABSTRACT 

Timing analysis of the INTEGRAL-IBIS and Swift-BAT light curves of the Su- 
pergiant Fast X-ray Transient (SFXT) IGR J16465-4507 has identified a period of 
30.32±0.02days which we interpret as the orbital period of the binary system. In ad- 
dition 11 outbursts (9 of which are previously unpublished) have been found between 
MJD 52652 to MJD 54764, all of which occur close to the region of the orbit we regard 
as periastron. From the reported flux outbursts, we found a dynamical range in the 
interval '^30-80. Although in this regard IGR J16465-4507 cannot be considered a 
classical SFXT for which typical dinamical ranges are >100, still our reported values 
are significantly greater than that of classical persistent variable supergiant HMXBs 
(< 20), supporting the idea that IGRJ16465-4507 is an intermediate SFXT system, 
much like few other similar cases reported in the literature. 

Key words: Gamma-rays: observations - X-rays: binaries - X-rays: individual: IGR 
J16465-4507 



1 INTRODUCTION 

During the last few years, the IBIS instrument 
l|Ubertini ct al."20 0J) on board the INTEGRAL gamma-ray 
satellite (Winkler et al.l 120031 ) has played a key role in 
discovering a previously unrecognized class of high mass 
X-ray binaries which display a very unusual fast X-ray 
transient behaviour and host a massive OB supergiant 
star as companion: the so-cal led Supergian t Fast X-ray 

Sgucr a et al.l 120051 . 



Tran s ients (SFXTs; 



,2005; 



lin't Zandl , _„ .__ 

I2OO6I ; iNegueruela et al.l |200^ To date approximately 10 
show clear characteristics of being a SFXT (plus the same 
number of candidates), characterized by long periods of 
low X-ray activity level (with luminosities or upper limits 
typically in the range lO'^^-lO'^^ erg s~^), interspersed with 
short, strong, flares lasting from a few hours to no more 
than a few days and reaching peak luminosities of 10^^-10^^ 
erg s~^. The identification of periodicity in the SFXT 
outburst behaviour, most likely representing the orbital pe- 
riod, is a key diagnostic for studying the geometry of these 
systems and hence for understanding the still unknown 
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phy sical mechan isms behind their very unusual behaviour 
(see ISidolil 12003 for a review). Because the fast outbursts 
from SFXTs are relatively rare occurrences, the discovery 
of such periodicities is a very challenging task which 
has been fulfilled only for an handful of SFXTs, mainly 
thro ugh long-term temporal mon i toring by INTEGRAL 
(e.g. IClark et al.' '2009'; Bird ct al. 2009; Jain ct al. '2009|; 
Zurita Heras et al. 2009 ; Sguera et al. 2007 ; Sidoli et al.l 
2007^ . 



IGR J16465-4507 classification as an SFXT is debat- 
able. Its initial discovery on 6 September 2004 was triggered 
by the INTEGRAL observation in the Norma region of a 
faint ~9 mCrab source (18-60 keV) which produced a strong 
~28 mCrab flare on 7 Septe mber 2004 before return ing to its 
undetectable pre-flare level (JLutovinov et al.ll2004h . demon- 
strating its fast X-ray transient nature. However since then, 
no more outburst activity has been reported from the source 
and its temporal and spectral behaviour above 20keV is 
largely unknown. A XMM-Newton observation performed 
a few days after the source discovery allowed the detec- 
tion of high intrinsic absorption (~6xl0^^ cm~'^) as well 
as a ~228s period wh ich was interpreted as the spin pe- 
riod of a neutron star ([Walter et al.l 120061 : iLutovinov et al.l 
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Figure 1. IBIS/ISGRI long term light curve (18-60 keV) of IGR 
J16465-4507. Each data point represents the average flux dur- 
ing one ScW (~2000s). The six activity periods listed in table [T] 
are indicated in the light curve by the numbers from 1 to 6. A 
zoomed view of the activity periods 2 and 3 (from MJD 53225.73 
to 53259.48) is shown in the inset. 



[2OO5I). A single bright infrared object falls in the refined 
XMM-Newton position (Zurita Heras & Walter 2004), opti- 
cal and infrared spectroscopy classified it as a supergiant 
star located at a distance of 9-5'^l\^ kpc, alth ough different 
spect ral types h ave been proposed (Q9 . 5Ia bvlNespoh et al.l 
20081 or B0.5I bv'Negueruel a et al.ll2006l : iRahoui et al.ll2008l ). 



Walter fc Zur ita Heras (20o3) reported a 22-50 keV quies 



cent flux of 7.7 x 10 erg cm s and suggest that IGR 
J16465-4507 is not a c lassical SFX T (for which the dy- 
namical range is >100, IChatvl l2010r ) but a kind of inter- 
mediate SFXT still having a dynamic range greater than 
that of a classical persistent variable supergiant HMXB 

i for which the dynamic r ange is <20) such as Vela X-1 
Krevkenbohm et al.ll2008ll. In comparison a very recent pa- 
per by iLa Parola et al.l ((20101) find no evidence of flaring 
activity from IGR J16465-4507 and inferred a dynamical 
range <10. 

Here we report a comprehensive temporal study of 
IGR J16465-4507, based on long term INTEGRAL moni- 
toring, which led to the discovery of 9 previously unpub- 
lis hed outbursts as well a s the ~30 days periodicity found 
bv lLa Parola et al.l (|2010l ). which we interpret as the binary 
period. 



2 DATA ANALYSIS 

Archival IBIS data from MJD 52652 to MJD 54764, 
giving an on source exposure of ~ 7.4 Ms, were pro- 
cessed with the INTEG RAL Off-hne Science Analysis (OSA 
iGoldwurm et al.l |2003| ) version 7.0. A hght curve for IGR 
J16465-4507 was generated on science window timescales 
(~ 2000 s) over the 18-60 keV energy range (flg. [ij. The 
light curve was then searched iteratively for outbursts. The 
shortest timescale region of the light curve that produces 
the highest detection signiflcance greater than 4 cr is found. 
This is the brightest outburst in the light curve. This section 
of the light curve is then removed from the total light curve 
and the next highest signiflcant section is searched for. 




Figure 2. Phase Dispersion Minimization of the outbursts found 
shows a minimum at 29.88 days with a 23% duty cycle. 




Figure 3. Periodogram for the IBIS light curve showing a period 
detection at 30.32 days. The side lobes to the period can be shown 
to be due to the window function of the light curve, while the 
lower significance peak at ~0.07 days~^ is an alias of the main 
peak. 



In total 11 outbursts have been found in the light curve 
(Table [1}. The outbursts range from short flares of 6.8 h to 
periods of flaring activity up to 87 hours. However, many 
of these outbursts are close together and if considered to be 
part of an overall period of flaring activity, only 6 periods 
of flaring remain. By folding the outburst times on differ- 
ent periods, an estimation of the duty cycle required over 
the orbit to explain all the outbursts can be made. This 
Phase Dispersion Minimization (PDM) technique produces 
a minimum of ~23% duty cycle at a period of 29.88 days 
in our data (Fig. [2]), which could indicate the orbital period 
and that outbursts occur over ~23% of the orbital phase. 
A Monte-Carlo simulation of the PDM performed by ran- 
domising the outburst time used 200,000 times within the 
length of the outburst provides an error estimate of ±0.46 d. 
This is in aggreement with the recent period found in the 
Swift-B AT light curve l|La Parola et al.|]2010l ). 
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Table 1. Outbursts found in the 18-60 keV ISGRI light curve, ordered by MJD. Peak luminosities have been calculated by assuming a 
distance of 9.5 kpc. It is clear due to the separation between the outbursts found that several could be combined into periods of flaring 
activity, indicated in column 1. 



Activity 
Period 



Sig 



Start 
MJD 



Stop 
MJD 



Length 
Days 



Peak Flux 



10 



-11 



erg cm 



Peak Luminosity 
lO^^ergs^i 



Orbital 
Phase 



1 


4.05 


53081.12 


53081.54 


0.40 


42.12 


2 


5.27 


53225.73 


53226.56 


0.83 


33.27 


2 


14.03 


53226.59 


53227.16 


0.54 


39.31 


2 


4.69 


53227.88 


53229.90 


1.99 


20.49 


2 


7.71 


53231.89 


53232.27 


0.36 


35.10 


3 


16.97 


53255.24 


53255.67 


0.42 


63.18 


3 


5.37 


53259.19 


53259.48 


0.28 


24.36 


4 


4.01 


53468.18 


53471.82 


3.63 


25.41 


5 


5.35 


53617.27 


53618.64 


1.34 


26.18 


5 


4.47 


53620.30 


53620.91 


0.59 


23.93 


6 


4.43 


53648.85 


53650.54 


1.68 


25.27 



4.52 


0.24 


3.57 


0.02 


4.22 


0.04 


2.20 


0.11 


3.77 


0.21 


6.79 


0.98 


2.62 


0.11 


2.73 


0.06 


2.81 


0.93 


2.57 


0.02 


2.71 


0.98 




Figure 4. IBIS 18-60 keV light curve folded on the discovered 
30.32 day period with the exposure for each bin overlaid (dashed). 
An ephemeris of MJD 53736.2 has been determined to place the 
peak of the emission at phase 0.0. 
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Figure 5. Folded IBIS light curve (solid line) with outbursts re- 
moved (dashed line). Arrows indicate the phases at which follow 
up X-ray observations have been made by different X-ray instru- 
ments. 
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Figure 6. Periodogram for Swift orbit light curve showing a sig- 
nificant period at 30.23 days. 

3 PERIODICITY AND SPECTRAL ANALYSIS 

The science window light curve (containing 2911 indepen- 
dent data points) was searched for signs of periodicity using 
the Lomb-S cargle periodogram method (|Lomblll97a : IScarglel 
Il982l . ll989r ). The Lomb-Scargle periodogram is shown in Fig- 
ure [3] with a clear peak evident at a frequency of 0.033 
days~^, corresponding to a period of 30.32 days, agreeing 
with the value determined by the PDM. The significance 
of th e peak was con firmed by performing a randomization 
test (jHill et al.ll2005l ): the fiux points of the fight curve are 
randomly re-ordered and a new Lomb-Scargle periodogram 
generated, the distribution of the power of the largest peak 
is an indicator of the significance of detection. From 200,000 
simulations we estimate the 99.9%, 99.99% and the 99.999% 
significance levels which are shown in Figure|3] note that the 
maximum power achieved in any of the randomized fight 
curves was <22. 

The error on the identified period was estimated using 
a Monte-Carlo simulation. Each flux measurement was ad- 
justed using Gaussian statistics within its individual error 
estimate to generate a simulated fight curve of the source. 
The corresponding periodogram was produced and the lo- 
cation of the maximum peak near the frequency of 0.033 
days~^ was recorded. From 200,000 simulations we estimate 
the period and its error to be 30.32±0.02 (la) days. 
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Figure |4] shows the phase folded hght curve using a pe- 
riod of 30.32 days and using an ephemeris of MJD 53736.2 
to centre the fold at phase 0.0. Assuming any enhancement 
to the emission to coincide with the neutron stars closest 
approach to the donor, we consider phase 0.0 to be perias- 
tron. To see the effect of the outbursts on the phase folded 
curve a new fold was performed excluding all of the out- 
bursts from Table [l] This second phase folded light curve 
is shown in Figure [S] and shows no substantial difference 
in shape, demonstrating that there is an underlying period- 
icity in the flux aside from the outbursts when the entire 
light curve is considered. The public light curves for IGR 
J16465-4507 from the Swift-3AT mission were also ana- 
lyzed using the Lomb-Scargle periodogram with the same 
tests used on the INTEGRAL light curve to see if there was 
any indication of the 30.32 day periodicity. The Swift orbit 
light curve showed a period of 30.23±0.04 d with a power 
of 22.43 (Fig. O, compatible with the period found in the 
INTEGRAL light curve. 

There is some evidence for an eclipse (Fig. |4]) in the 
phase folded INTEGRAL light curve. However, this cannot 
be completely disentangled from low exposure in the relevant 
phase bins and the statistical significance of this dip is low. 
The feature is also not apparent in the Swift phase folded 
light curve leading to the conclusion that this is an exposure 
artifact and not a real feature. 

The two outbursts with the highest significance in Ta- 
ble [T] (2 and 3) are shown in the inset of Figure [T] Such 
IBIS data have been unsuccessfully searched for pulsa- 
tions. The brightest flares cannot be seen by JEM-X ei- 
ther on science window timescales or in deep mosaics in 
the 4-20 keV energy range, where the science windows ex- 
posure of 2.8 ks was not sufficiently long to detect the 
source. This results in an upper limit for the fiux dur- 
ing the outburst of 8x10^^^ ergcm^"^ s^^. Only one Swift- 
XRT observation is available with only 69 events in 2ks 
and so it has not been possible to extract a meaningful 
spectrum and to search for pulsations. However this does 
give an estimated 0.2-10 keV absorbed (unabsorbed) fiux 
of 9.1x10"^^ erg cm"^s"^ (4.9xl0~" ergcm"^ s"^), based 
on a power law index of 1.0, a Galactic hh ~ 2.12x10^^ 
and an intrinsic uh = 6x10'^^. The flux has also been mea- 
sured by XMM-Newto n as 3±lxl0~^^ ergcm"^ s~ ^ in the 
4-10 keV energy range (|Zurita Heras fc WalterlbOO^ I and as 
8.95t3;i X 10"^^ ergc m~^s~^ in the 0.2-10 keV energy range 
using Suzaku-XIS (jMorris et al.ll2009y These measurements 
approximately follow the drop in flux in the phase folded 
light curve (Fig. [S]). 

Spectra for the two outbursts with the highest sig- 
nificance detection (n. 2 and 3 in table 1) have been 
extracted from the IBIS data. The 18-60 keV spectrum 
of the discovery outburst (n.3) is reasonably fit by a 
bremsstrahlung model (x?=0.6, 15 d.o.f., kT=18.6l^ 5 keV, 
flux equal to 1.2x10"^" ergcm"^ s~^) or alternatively by 
a power law (x2=0.7, 15 d.o.f, F = 3.0±0.3). Similarly, 
the newly discovered outburst (n.2) is also reasonably fit 
by a bremsstrahlung (x^=l-2, 15 d.o.f., kT=29t" keV, 
flux equal to 9.5x10"^^^ ergcm"^ s~^) or by a power law 
(X?=1.24, 15 d.o.f.,F = 2.5±0.3). It is evident that both 
outbursts have a similar soft spectrum and a comparable 
average fiux of activity, furthermore the reported spectral 



shapes are very similar t o those seen during outb urst activ- 
ity of other SFXTs fe.g. lSguera et al.ll2008l . l2007l ). 

Taking into account the 18-60 keV bremsstrahlung 
spectral shape of both previous outbursts, we extrapo- 
lated their 0.2-10 keV soft X-ray fiux obtaining a value 
of ~3xl0~^*' ergcm"^ s~^. To date the lowest and high- 
est measured unabsorbed soft X-ray fiuxes are equa l 
to ~3.6xl0"^^ergcm"^s"^ (|Zurita Heras fc Wahej |2004| '1 
and ~4.9xl0~^^ ergcm"^ s~^ (this paper) so the source dis- 
plays a dynamical range from 6 to 83. As for the hard X- 
ray band, if we compare the 22-50 keV quiescent fiux (~ 
7.7 X 10~^^ergcm-^s- MWalter fc Zurita Herad l|2007l )) to 
the outbursts' fluxes as listed in table 1, then the dynam- 
ical range is from 27 to 82. We note that such values are 
significantly higher th an that recently reported (<10) by 
iLa Parola et al.l (|20ld ^ as based on 54 months BAT mon- 
itoring (15-50 keV) during which no evidence for fiaring 
activity was found. However, we point out that such appar- 
ent dynamical range incongruences can be explained by the 
different BAT and IBIS instrumental capabilities. IBIS is 
particularly suited to detect short outbursts and so to re- 
veal the dynamic range of fast X-ray transient sources, such 
as IGR J16465-4507, thanks to its higher instantaneous sen- 
sitivity for shorter observation lengths (i.e. science windows, 
~ 2000 s duration) which match very well the duration of the 
fiares. BAT is more suited to long-term monitoring thanks to 
its more continuous coverage but could eventually miss the 
detection of short fiares which, if close to the IBIS sensitiv- 
ity, are very likely not detectable by BAT due to its poorer 
instantaneous sensitivity. It is worth noting that IBIS and 
BAT show to be complementary and their combination is 
important in investigating the properties of SFXTs. 

In order to test the recurrenc e of emission at p eriastron 
we carry out a recurrence analysis (jBird et al.l2009l ). For this 
test we calculate the light curve significance for a region of 
4 days either side of the points considered as periastron and 
apastron (approximately 25% of the orbit). This will high- 
light if we see a low level emission at periastron passages 
compared to apastron. Comparing the distributions of sig- 
nificances for apastron and periastron using a Kolmogorov- 
Smirnov test we find a probability of 74% that they have 
been drawn from different distributions. 

4 DISCUSSION 

The orbital characteristics and stellar properties for IGR 
J16465-4507 are ve ry similar to those of SAX J1818. 6-1703 
(Bird et al. 2009; Zurita Heras et al. I [20091. The peak ffux 
from the phase folded light curve with the o utburst re- 
moved of SAX J1818. 6-1703 as determined by iBird et all 
(|2009l ) is ~2.2cs"^ in the ISGRI 18-60 keV band. Com- 
paring this to the peak flux of ~0.55cs~^ from the same 
light curve for IGR J16465-4507 shows that IGR J16465- 
4507 is ~4 times weaker if we assume that these systems are 
very similar. The difference in flux is compatible with the 
known distance measurements for both sources. However, 
SAX J1818. 6-1703 shows much less absorption with Nh = 
(6.0 ± 0.7) X 10^^ cm"^ compared to iVi/ ~ 6 x 10^^ cm"^ 
for IGR J16465-4507. 

Since masses and periods are known from the spec- 
tral type (|Drilling fc Landoltll2000l ) and assuming that the 
compact object is a neutron star of mass 1.4Mq, Kepler's 
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Figure 7. Corbet Diagram IjCorbet et al.ll2006r i siiowing tiie lo- 
cation of IGR J16465-4507 (clear inverted triangle) compared 
with other SFXTs (purple inverted triangles). Be X-ray transient 
systems (triangles), Wind accretors (boxes) and Roche-lobe ac- 
cretors (circles) are shown for reference. 



third law gives the semi-major axis for the system of 121.9- 
126.3-Rq. As we see no evidence for Roche-lobe overflow in 
the system the eccentricity of the system must be <0.8 for 
the smallest possible donor star size and <0.6 for the largest 
donor star size of 26Rq . 

It is important to point out that whereas classical 
SFXTs have a remarkable dynamical range (>100), we found 
for IGRJ 16465-4507 significantly lower values (in the range 
~30-80). Although in this regard IGR J 16465-4507 cannot 
be considered a classical SFXT, still its dynamical range 
deserves some attention since it is greater than that of 
classical persistent variable supergiant HMXBs (<20). Our 
findings on the source temporal behaviour and its dinam- 
ical range are in agreement with the results reported by 
Walter & Zurita Heras (2007) and strongly support their 
idea that IGRJ16465-4507 is an intermediate SFXT system, 
much like few other similar c ases reported in the literature 
l| Walter fc Zurita Herail2007l ). 

As this system has both the orbital period and spin pe- 
riod known then it may be plotted on the Corbet diagram 
(Fig.[7I). This system resides between the regions dominated 
by wind accreting systems an d Be X-ray transien t systems, 
much like IGR J18483-0311 jSguera et al.1 12007| ) and SAX 
J1818. 6-1703. When looked at in the context of the other 
SFXTs, IGR J16465-4507 sits in the middle of the class. 
This implies a longer orbit than typical wind accreting su- 
pergiant systems, but without evidence for accretion from 
a circumstellar disk, as in Be X-ray binaries. Our measure- 
ment of the orbital period in this system, combined with 
bursting behaviour of the source and constraints on the ec- 
centricity of the orbits support the picture of accretion from 
a clumpy wind (jNegueruela et al.l 120061 : iDucci et al.l 12009^ . 
The underlying modulation in the phase folded light curve 
can be explained by the eccentricity of the system as with 
IGR J17544-2619 (ciark et al.ll2009i ). However, as shown in 
figure[Sl an underlying periodic emission, around periastron, 
does appear to be present in the phase folded light curve 
even when bright outbursts have been subtracted. 



5 CONCLUSIONS 

IGR J16465-4507 has been shown to exhibit an ~30 day 
periodicity in both the INTEGRAL-lBlS and Swift-BAT 
light curves. This period is attributed to the orbital pe- 
riod of the system and corresponds to a semi-major axis of 
121. 9-126. 37?0. This system is remarkably similar to SAX 
J1818. 6-1703, except that IGR J16465-4507 is known to be 
further away, helping explain the lack of underlying emis- 
sion that can only be seen in the phase-folded light curve of 
the entire INTEGRAL data set and not in the recurrence 
analysis. A low dynamic range is observed, which is higher 
than that of a typical supergiant HMXB. This would sug- 
gest that IGR J16465-4507 is an intermediate SFXT and 
that these may form a continuum ranging from typical su- 
pergiant HMXBs to the extreme SFXTs. 
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